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Ah. St  r.ic t 

This  report  con.siders  an  array  of  simultaneous  under¬ 
water  nuclear  explosions  general  ini’  plane  tidal  waves  on 
the  continental  shelf.  Computer  simulations  are  used  to 
study  this  wave  generation  process.  They  are  based  on  non¬ 
linear  shallow  water  theory  with  an  initial  condition  of  a 
stationary  raised  cone  of  water.  'I'lie  sensitivity  of  the 
wave  lieight  to  different  yields,  combination  of  different 
yields,  spacing,  alignment,  timing,  and  bottom  slope  is 
determined.  The  reference  values  used  for  this  study  are  a 
yield  of  five  kilotons  TNT,  one  kilometer  spacing,  water 
depth  of  100  meters,  and  no  bottom  slope.  In  addition  to 
the  sensitivities,  the  amount  of  non-linear  effect  is  illus¬ 
trated  by  comparison  of  wave  contour  printer-plots  generated 
by  linear  and  non-linear  shallow  water  theory.  The  computer 
results  indicate  the  wave  generation  is  relatively  insensitive 
to  realistic  variations  in  the  parameters. 
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IM.ANi;  TIDAL  IvAVL.S  (’.liMI-RATTA)  BY 
AN'  ARRAY  Oi' 

SIMULTANLOUS  UNUi.RWATi:R  LXPLOSIONS 

I .  T  n t  roduc  t i on 


Background 

This  thesis  originated  from  the  Al'Ii  iLNIi-SOM  class 
design  project  which  analyzed  a  shallow  underwater  mobile 
CSUM)  basing  alternative  to  the  land  based  MX  intercontinental 
ballistic  missile  weanon  system  (Ref  1).  The  SUM  system  would 
be  deployed  over  tlie  continental  shelf  of  the  United  States. 

This  large  area  of  deployment  with  the  location  of  the  SUMs 
unknown  would  require  a  large  number  of  warheads  to  defeat  the 
weapon  system.  However,  an  easy  defeat  of  the  SUM  based  weapon 
system  lias  been  proposed  which  uses  a  plane  tidal  wave  generated 
by  an  array  of  simultaneous  underwater  explosions.  The  tidal 
wave  would  travel  across  the  continental  shelf,  disabling  the 
SUMs.  Computer  simulations  based  on  shallow  water  theory  verify 
the  formation  of  a  plane  tidal  wave.  The  wave  energy  resulting 
from  the  underwater  explosion  was  approximated  with  an  initial 
condition  of  a  stationary  raised  cone  of  water.  Further  study 
of  this  plane  tidal  wave  generation  is  tlie  topic  of  this  thesis. 

Problem  and  Scope 

Tlie  primarx"  puriiose  oP  this  stiuly  is  to  determine  the 
sensitivity  of  tlie  iflanc  tidal  wave  generation  to  the  explosions' 
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yield,  yield  coiub  i  nat  ioi\s  ,  spaciai\  between  explosions,  align¬ 
ment  of  explosions,  timing  of  explosions,  and  slope  of  the 
continental  shelf.  The  reference  values  used  are  five  kilotons 
TNT  for  yield,  one  kilometer  spacing  between  explosions,  and 
a  water  depth  of  1(10  meters.  flie  yield  is  varied  from  2  to  8 
kilotons.  The  combination  of  two  yields  is  varied  as  a  ratio 
of  one  yield  to  the  other  yield  from  .25  to  1  with  an  average 
yield  of  5  kilotons.  The  spacing  is  varied  ±401.  from  the 
reference  value  of  one  kilometei’.  The  misalignment  is  varied 
from  0.0  to  0.5  kilometer  and  tliie  timing  is  varied  from  0.0  to 
5  seconds.  The  bottom  slope  is  varied  from  0,0  to  0,002, 

In  addition  to  the  sensitivities,  the  amount  that  the  non¬ 
linear  effects  contribute  to  tlie  generation  of  the  |)lane  tidal 
wave  is  illustrated, 

Assumpt  ions 

There  are  two  main  assumptions.  The  first  is  that  an 
explosion  on  the  ocean  floor  wliich  produces  a  large  enough  gas 
bubble  for  blowout  will  generate  sl\allow  water  waves.  The 
second  assumption  is  part  of  shallow  water  tlieory.  It  assumes 
the  water's  vertical  acceleration  doesn't  affect  pressure  and 
tlie  water's  vertical  velocity  iloesn't  affect  the  conservation 
of  momentum  equation. 

Approach  and  Set|uence 

The  general  approacli  starts  by  analysing  the  problem 


physically  to  determine  tlie  type  of  wave  generated 


realistic 


parameter  reference  values,  aiul  model  in;^  of  initial  condition. 
Next,  the  differencing  scheme  is  develojied  !)>■  deriving  tlie 
shallow  water  equations  and  differencing  the  equations.  The 
scheme  is  used  to  simulate  the  wave  generation.  The  wave 
heiglits  are  evaluated  to  determine  sensitivity  to  several  of 
the  parameters.  Finally,  the  non-linear  effects  are  illus¬ 
trated  . 


I  I  .  Ana  lysis  o  1  I  *  r  o  h  1  c'lii 

Tlie  analysis  of  this  piohlcm  focuses  on  tlirce  areas. 

First  is  tlie  rationale  and  basis  for  usini>  shallow  water 
theory.  Next  is  tlie  determination  of  the  reference  values 
for  the  indeiiendent  variables  and  tlieir  rant>e  of  variation. 

■fhe  last  area  is  tlie  initial  condition  and  its  relationship 
to  the  explosions'  yield. 

Rat  io n-  le  for  Shallow  Water  The o r >' 

The  rationale  for  usiny  sliallow  water  theory  is  primarily 
based  on  the  type  of  cavity  formed.  Shallow  water  theory  can 
be  used  if  the  total  depth  of  tiie  water  moves  horizontally 
together.  The  explosions  used  for  this  problem  generate  a 
gas  bubble  witli  a  radius  greater  than  the  deptli  of  water. 

Since  tlie  expanding  butible  will  horizontally  excite  the  total 
depth  of  tlie  water,  shallow  water  waves  arc  expected.  On  this 
basis,  shallow  water  theory  is  used. 

Research  at  l.os  Alamos  Scientific  Laboratory  (LAST)  found 
that  shallow  water  theory  was  not  applicable  for  a  near  surface 
explosion  fRef  (i).  However,  the  difference  in  burst  location 
can  account  for  this.  The  near  surface  burst  used  by  LAST 
formed  a  hemispherical  surface  cavity  with  a  radius  equal  to 
l/b  the  depth  of  the  water.  The  surface  cavity  would  not  hori¬ 
zontally  excite  the  total  depth  of  the  water  at  the  sanie  time. 
Therefore,  shallow  water  waves  are  not  expected,  which  is 
consistent  with  shallow  water  theory  not  accurately  duplicating 
the  experimental  data. 
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Parameters  Re f e  reaee  Vaiues 

The  reference  values  of  ilte  i  ml  epeiuk  n  t  pa  rajiiet  e  r  s  and 
rani>c  of  variation  were  selected  based  on  many  factors.  The 
water  depth  and  bottom  slope  are  based  on  the  averat^e  depth 
and  slope  of  the  continental  slielf.  The  yield  was  selected 
based  on  wave  generation  efficiency  and  the  minimum  yield 
to  ensure  blowout.  The  spacing  was  limited  by  the  accuracy 
of  the  bomb  [placement.  'I'he  alignment  and  timing  are  assumed 
perfect  for  the  leference  value. 

Depth  and  Slope.  The  reference  water  dejitli  is  al)out  the 
average  of  the  continental  shelf's  depth  (0  -  bOO  ft).  i  lie 
water  deptli  reference  value  is  100  meters  witli  no  variation 
except  when  bottom  slope  is  used.  The  average  bottom  slo))e 
of  tlie  continental  shelf  is  .0002  to  .  0004  based  on  a  width 
of  300  to  600  miles.  The  reference  value  for  bottom  slope  is 
zero  with  a  variation  of  0.0  to  .002. 

i'iming  and  Alignment.  The  reference  value  for  timing  is 
zero  seconds  with  a  range  of  variation  of  0  to  .3  seconds. 

The  reference  value  for  alignment  is  also  zero  with  a  range 
of  variation  of  0  to  500  meters.  Misalignment  distance  is 
the  perpendicular  distance  from  tlie  array  centerline  to  tlie 
misaligned  explosion. 

Spacing.  The  spacing  reference  value  is  limited  at  the 
lower  end  by  the  accuracy  of  tlie  warhead  placement.  According 
to  S  c  L  c  n  t i f i c  Amor i can  ,  the  1085  estimated  circular  error 


pronahli'  ((lliPI  (nr  '’iissian  missiles  is  .M  kilometer  (Ref  2:541. 
The  rofereiu'c  spaciiis;  was  limited  to  a  distance  of  4*(d;P. 

Based  on  tliis,  tlie  reference  v'alue  of  one  kilometer  is  used 
w  i  t  li  a  va  1'  i  a  t  i  on  of  .  (>  to  1.4  k  i  1  oiiie  ter. 

Yield.  file  reference  yield  (Y)  is  based  on  the  efficiency 
of  wave  formation  and  tlie  minimum  yield  required  for  a  bloi\Jout 
case.  I'he  efficiency  as  a  function  of  yield  is  indicated  by 
the  wave  heiyht.  I'he  wave  heiylu  (li)  is  proj^ort  ional  to 
(Ref  5:133).  I'lierefore,  tiie  lowest  possilile  yield  will  be  used 
to  obtain  the  hii>liest  efficiency. 

The  minimum  yield  required  to  produce  a  maximum  hemisper- 
ical  bubble  with  a  radius  equal  to  tlie  water  depth  of  100  meters 
is  1 .  kilotons.  This  is  based  on  the  following  ajiprox  imat  ions 
and  calculations.  The  energ>'  required  to  form  tlie  hem  i  splier  i  ca  1 
bubble  is  approximated  by  the  energy  rcuiuired  to  form  an  equal 
volume  spherical  bubiile  with  coinciding  centroids.  Approximately 
half  of  the  explosion's  energy  is  absorbed  by  the  ocean  floor 
(Ref  5:126).  Using  these  a|iprox  imat  ions  and  IA|  (1)  for  the  maxi¬ 
mum  bubble  radius  (Ref  5:113),  tlie  minimum  yield  is  determined 
as  follows.  I.  is  the  max  bubble  radius,  Y  is  tlie  yield,  1) 
is  tile  depth  of  burst,  R^  is  the  radius  of  the  spherical 
bubble,  and  Rj^  is  the  radius  of  tlie  hem  i  sjihe  r  i  c  a  1  bubble. 
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70.4  meters 
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To  obtain  the  minimum  yielJ,  use  l;q  (3)  with  L  =  , 

D  equal  to  the  depth  of  the  hemisphere's  centroid,  and 
Y  =  ^  Y  to  account  for  the  energy  absorbed  by  the  bottom. 

j  Y(kt)  =  (62.5  +  10) 

Y  =  1.7  ktons 


Tlicrefore,  the  yield  should  be  equal  to  or  greater  than  tlie 
minimum  yield  of  1.7  kilotons,  A  reference  yield  of  5  kilotons 
was  selected  since  it  is  small  and  vet  more  than  twice  the 
minimum  \'ield  to  ensure  shallow  water  waves. 

Initial  Condition 

The  initial  condition  used  to  represent  the  wave  energy 
imparted  by  the  explosion  is  very  simplistic.  It  consists  of 
a  raised  stationary  cone  of  water.  I'he  cone  surface  has  a 


constant  slope  of  0.2  with  its  volume  JepeniJent  on  the 
yield. 

The  relationship  between  t  lie  yield  and  tlie  potential 

energy  of  the  cone  is  based  on  the  scaling  laws  for  the 

blowout  case  (Ref  5  : 1 32  - 1  .  In  the  blowout  case,  the 

wave  height  is  scaled  by  .  Since  the  t>ot  t  i i  energy 

of  the  wave  is  proportional  to  the  wave  heigiit  squared,  it 

is  scaled  by  y^/^  ,  Therefore,  the  potential  energy  (PT:) 

1/3 

of  the  cone  is  equal  to  some  constant  times  Y 

Pii  =  Const*Y^'^^  (4J 

c  on  e 

The  constant  is  obtained  liy  analyzing  tlie  potential 

energy  of  the  cavity.  Starting  with  the  minimum  yield  of 

1.7  kilotons,  the  maximum  bubble  radius  (1.)  is  100  meters. 

Assume  the  water  displaced  by  the  bubble  is  now  located  at 

the  surface.  The  amount  of  increase  in  potential  energy  is 

determined  as  follows;  where  p  is  water  density  (1  000  kg/m"^i 

2 

and  g  is  the  acceleration  of  gravity  (9.8  m/s''}. 

APP  =  (62.5)(jJti  L^)  pg  =  1.28*10^^  Joules  (5) 

As  the  bubble  collapses,  much  of  the  potential  energy  is 
lost  in  turbulence.  One  half  of  the  wave  energy  is  assumed 
to  be  lost  in  turbulence.  Based  on  this  assumption,  the  constant 
is  determined  as  follows; 
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PO  (J)  =  5.4*10^^  Y^/^(kt) 

c  on  e 


(7) 


The  size  of  the  cone  is  determined  by  the  amount  of  wave 
energy  imparted  by  the  explosion.  The  potential  energy  of 
the  cone  of  ivater  is  set  equal  to  the  energy  obtained  from 
Eq  (7).  Solving  for  the  height  (h)  of  the  cone's  peak 
determines  the  cone's  size  with  a  constant  slope  (s)  of  0.2. 
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III. 


Clomput  or _ Prop,  ram  Dev e  1  opiiiLii t 

The  computer  propram  ilevel  ofHiicn  t  includes  the  derivations 
and  differencing  of  the  equations  and  treatment  of  tiie  boundary 
conditions.  Both  the  non-linear  and  linear  equations  will  be 
derived  and  differenced.  In  addition,  a  check  of  the  differ¬ 
encing  scheme  is  accomplished. 

Derivation  of  r.quations 

The  sliallow  water  wave  equations  are  derived  from  the 
conservation  of  mass  and  conservation  of  momentum  equations 
using  the  3 -d imens iona 1  rectangular  coordinate  system.  Die 
viscosity  and  compressibility  of  the  water  arc  ignored.  'Die 
coordinate  system  is  as  shown  in  Figure  1  with  the  free-standing 
water  surface  at  z  =  0  .  The  actual  water  surface  position 

is  described  by  ll(x,y,t)  and  tlie  bottom  is  descrilicd  by 
b(x,y)  .  The  averaged  velocity  components  u  ,  v  ,  and  w 
are  in  the  x  ,  y  ,  and  z  directions,  respectively,  witli  S 
representing  the  velocity  vector  (u,v,w) 

The  derivation  starts  with  the  conservation  of  mass  equation 
(Ref  4:2,24)  and  its  integration  witli  respect  to  z 


Dt 


p  V  •  -S 
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(in 


9u  ^  3v  ^  9w  ^  „ 

9x  9y  9z 


(12) 
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Now  the  first  approximation  must  he  made  for  the  shallow 
water  theory.  The  vertical  acceleration  of  the  water  is 
assumed  not  to  affect  the  j)ressurc  (P)  (Ref  8  :  23-24).  This 
results  in  hydrostatic  pressure,  P  =  updl-z)  where 
b  -  z  ^  H  .  Hydrostatic  pressure  provides  a  horizontal 
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pressure  gradient  whieli  is  i  tulcpendeiit  of  z  .  Therefore, 
the  horizontal  aeeelerat  ioys  and  velueities  are  independent 


of  2  .  With  tliis  ajjprox  i  mat  i  on  ,  t  lie  integration  of  the 

conservation  of  mass  equation  (17)  is  completed. 
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(19) 


The  conservation  of  momentum  otiuation  (Kef  4:3,24)  is 
used  to  derive  the  velocity  equations.  The  prev'ious  approxi¬ 
mation  is  expanded  to  accomplish  this,  Tlie  vertical  velocity 
w  is  assumed  to  have  a  negligible  effect  on  the  conservation 
of  momentum.  Therefore,  the  velocity  vector  S  becomes  (u,v)  . 
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The  ec|uation  for  v  is  obtained  similarly. 
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Hquatiorib  (19),  (21)  and  (22)  are  the  non-linear 
shallow  water  equations.  liie  linear  shallow  water  equations 
are  obtained  by  assumini>  u  ,  v  ,  II  ,  and  their  derivatives 
are  small  enough  so  their  scjuares  and  products  can  be  neg¬ 
lected  (Ref  8:24).  This  results  in  linear  equations,  where 
-b  is  simply  the  depth  of  the  water. 
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(24) 


(25) 


If  tlie  water  depth  is  assumed  constant,  the  variables 
u  and  V  can  be  eliminated  as  follows  to  obtain  the  linear 
wave  hq  (30)  (Ref  8:25). 
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ny' 


This  shows  that  the  cliaiac  t  er  i  s  t  i  c  uavo  velocity  is  v'gC-bJ 


Uif Terence  Schei.ie 


The  second  order  explicit  diTTerence  scheme  is  three 
levels  in  time  and  space,  conservative,  and  lias  an  amplifi¬ 
cation  factor  of  unity.  Correct  space  and  time  centering 
provide  an  amplification  factor  of  unity  if  the  Courant- 
Fr iedrich-Lewy  criterion  is  satisfied  (Ref  7:288).  This 
criterion  is  -  *  •  Using  the  charac¬ 

teristic  velocity  equally  distributed  between  u  and  v  , 
the  criterion  value  is  0.44  for  At  =  . 5  sec  and  Ax  =  Ay 
=  50  meters  .  The  difference  equations  are  based  on  the 
integral  conservation  of  mass  ITi  (19)  and  tlie  conservation 
of  momentum  liqs  (21)  and  (22). 
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Kinematic  viscosity  can  be  added  to 
equations  to  account  for  viscosity  or  to 
instability.  I'lie  kinematic  viscosity  is 


the  last  two 
help  prevent 
incorporated  by 
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adding  to  tiu'  riglit  side  of  l.(is  (dJ)  and  (?ia)  a  kinematic 
viscosity  coefficient,  multiplied  In  tlie  l.aplacian  of  u 
and  V  ,  respectively.  A  slight  diffusion  coupling  might 
also  he  added  to  keej)  the  independent  conservative  flows 
in  step.  This  is  done  by  adding  llte  pioduct  of  a  coupling 
coefficient  and  tlie  Laplacian  of  11  to  the  right  side  of 
liil  (31). 

■file  difference  eciuations  use  the  indices  i  ,  j  , 
and  n  for  tlie  variables  x  ,  y  ,  and  t  ,  respectively. 
The  equations  are  second  order  in  time  and  space  with  no 
kinematic  viscosity  or  diffusion  coupling  as  shown  lie  low. 
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The  linear  d  i  f  Cc  fence  e(|iiation  is  based  on  tlie  linear 
wave  lii|  (30).  It  is  second  order  with  space  and  time 
centering. 
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The  difference  scheme  used  a  stantlard  rectangular  mesli 
with  all  the  point  values  calculated  for  each  time  step. 
Syiiiinetry  w'as  used  as  much  as  possilile  with  the  aid  of 


reflecting  boundaries. 


Buiuid a ry  doiul  it  ions 

Ke  I'l  ec  t  i  bouinlai)'  (.oiul  i  t  i  on  s  a  lo  usui  in  tiic  proi^rain 
to  take  advantage  of  symmetry  (Ref  7  :  d  iS.S  -  J  d  1  )  .  I',  i  t  h 

reflecting  boundaries,  tlie  boiuuiar)'  conditions  are  set  up 
to  ensiii'e  no  mass  floiv'  aouss  the  !>onndaiy. 

first,  tlie  velocity  components  perpendicular  to  tlie 
boundary  are  set  to  zero  at  the  boundai-y  and  botli  velocity 
components  are  zero  at  the  coi-ners.  1  he  perpendicular 
velocity  components  just  outsiile  the  boundary  are  set  epual 
to  the  negative  of  the  pe  I'jKMid  i  cu  1  a  r  components  Just  inside 
tlie  boundary.  The  parallel  velocity  coinjionents  just  outside 
the  boundai')'  arc  set  equal  to  the  parallel  velocity  component 
just  inside  the  boundary.  The  surface  height  just  outside 
the  boundary  is  set  equal  to  the  surface  height  just  inside 
the  boundary.  The  corners  are  treated  similarly  with  the 
corner  [joint  being  the  center  of  symmetry  (e.g.,  if  corner 

[joint  is  (2,2)  ,  the  V(l,2)  -  -V(5,2)  ,  11(1,2)  =  0  , 

and  U(1  ,  1)  =  -  U(a,5)  1  . 

V^'  i  f  i  cat  ion 

■fhe  [jrogram  based  on  these  difference  et[uations  and 
boundary  conditions  was  verified  by  du]Jlicating  the  results 

obtained  by  Mader  using  bis  SiVAN  code  (Ref  b)  .  The  problem 

solved  by  Mader  consisted  of  a  .5  meter  radius  hem  i  s]Jbe  r  i  ca  1 
surface  cavity  in  water  three  meteis  dee|j.  The  cavity  was 
the  result  of  a  near  suil'ace  burst.  This  [Jiogram  used  the 
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LompiUcr  s  1  iiai  1  a  t  u>ii  s  hasL-il  on  t  ho  il  i  1' 1  e  i  ciicc  ci|Uulions, 
initial  o  oiul  1 1  i  on  s  ,  aiul  houndary  coiulilions  wore  Lisod  to 
dotei'iii  ino  tlie  sensitivity  ol  the  tidal  ivave  i'.ene  r  a  t  i  on  . 
ligure  1  illustrates  tlie  yr-iicial  shape  o  i'  tlie  tiLlal  wave 
generated  h\  the  computer  simulations,  Tlie  wave  velocity 
is  33(m/s),  which  is  ver_\  close  to  the  clia  lac  te  r  i  st  i  c  v>ave 
velocity  of  31.3(m/s).  i'he  initial  coiulition  used  produced 
a  wave  with  a  period  of  about  50  seconds  wliicii  is  nearlv 
thi'ce  times  the  exj'ected  period  of  IH  secoiuis  based  on  tlie 
re  1  a  t  ion  sli  i  [1  I’eriod  ^  11.1  {'lef  ihis 

indicates  tlie  wavelength  of  IbOO  meters  produced  by  this 
initial  corulition  is  also  about  tliree  times  too  long.  Tlie 
height  of  the  initial  wave  aliove  the  surface  was  selected 
as  the  indicator  of  the  vvave  generation  sensitivity.  The 
wave  height  is  measured  at  a  nearly  constant  distance  of 
8  kilometers  from  the  explosions.  ITie  wave  height  was 
obtained  from  the  p  r  i  n  t  e  r  -  p  1  o  t  s  which  expresseil  the  wave 
height  within  increments  of  .35  meters.  The  error  bars  on 
the  graphs  indicate  the  variation  in  wave  height  along  the 
plane  tidal  wave.  both  the  linear  and  non-linear  difference 
schemes  used  the  same  time  step  of  0.5  seconds  and  cell 
size  of  50  iiu'ters  square. 


Figure  2.  Illustration  of  Plane  Tidal 


V  I  c  1 J  Sell  sit  1  V  i  t  y 

The  wave  lieij^^ht  sensitivity  to  yieltl  is  proportional 
to  .  This  does  not  agree  with  the  sealing  of  wave 

height  as  used  for  tlie  initial  condition  (Ref  5:133). 

However,  it  does  come  close  to  the  scaling  of  wave  lieight 
as  [oj.  titjL'p  water  and  for  shallow  water  used 

by  lilasstone  (Ref  3:272  273).  I'he  sensitivity  to  yield  over 
tlie  range  of  2  -  8  kilotons  is  +4 . 5  (m/ (  kt  )  ^'^'^)  . 

The  wave  sensitivity  to  combinations  of  two  yields 
whose  average  is  5  kilotons  is  small.  i'he  sensitivity  is 
determined  relative  to  tlie  ratio  of  the  small  yield  over  tlie 
larger  yield.  As  shown  in  figure  4,  tlie  wave  lieiglit  is  not 
very  sensitive  to  combinations  of  yields  over  a  practical 
range  of  0,5  to  1. 

Spa  c  Ml  g  S  e  n_s  1 1  j  v  ji  ty 

i'he  spacing  sensitivity  was  determined  over  tlie  range 
of  .6  to  1.4  kilometers.  The  wave  height  did  not  cliange 
linearly  with  spacing  as  shown  in  1 igure  5.  Hver  the  range 
of  0.8  -  1,2  Kilometers  s/wicing,  the  apjiroximate  sensitivity 
is  -11.2  5(m/km)  sii.icing. 

A  1  ygniuen  t  Sen  sit  i  v  1 1  y 

i'he  sensitivitN'  to  misalignment  w.!--  termined  for 
misalignments  up  to  .5  kilometev  a--  indicated  in  I'igure  h. 
The  misalignment  situation  w.is  set  up  b)-  h.iving  every  otlier 
explosion  offset  a  riven  dist.ince  from  the  array.  flic 
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rigure  5.  Yield  Sensitivity 
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rijiure  4.  Yield  Combination  Sensitivity 
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til  i  s;i  1  1  giiiiieii  t  Jist.iiiii.  i  ililiiicvl  .is  this  dIIscI  Jistancc. 

The  a  ppi  o\ inui  t  e  sen  ^  i  t  i  v  1 1  v  over  the  raii).',e  ol  0  to  .5 
kilometer  is  -J.7(ui/km)  m  i  sa  1  i  jpiineii  t  .  Tlic  m  i  sa  1  i iiiiient  s 
did  iiierease  tiie  wave  lieii'ht  v.iriation  .iloiii;  the  wave  front 
as  indicated  liy  t  lie  eiior  ha  i  s  . 

I  iiii  i  111,  Sens  i  l  i  v  i  t  y 

Tlie  wave  heii’lit  sensitivity  to  tnniiut  delays  o  f  0  -  5 
seconds  is  very  small,  as  shown  in  I  ipnre  7.  Tlie  time 
delay  situations  are  ino>.leled  hy  Iniviny  ever)  other  explosion 
delayed  the  same  amount  i.)  f  time.  In  the  ranyte  of  ('  -  d  • -S 
seconds  delay,  the  seasitivit)  is  only  - . Od h ( m/ sec )  time 
delay. 

ho  t  t  oiii  S 1 0|K)  Sen s  i  i  i  \  it) 

The  average  slope  of  tlie  continental  slielf  had  a 
negligible  effect  on  the  w.ive  lieight  in  8  kilometers  of 
travel,  fhe  liottom  was  treated  as  being  smootli  and  with 
constant  slope.  llie  slope  was  varied  from  0.0  to  O.OOd 
as  indicated  hy  I  igure  .S .  The  wave  heiglit  sensitivity  at 
8  kilometers  from  tlie  explosions  is  rl8d(m)  per  the  slope. 
Assuming  the  ell'ect  was  constant  over  the  8  kilometers, 
the  sensitivity  was  +d  .1  { m/ km )  per  kilometer  traveleil  by  the 
wave  per  the  slope.  Tor  example,  using  the  averaged  conti¬ 
nental  shelf  slojie  of  ().0l)0.>,  the  cxiiected  increase  in  wave 
height  in  100  kilometers  is  0.7  meters. 
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Noil  -Li  near  hi'  fee  t  s 

The  non-linear  efrects  arc  i  1  lust  rat  cil  l)y  eoiiipa  r  ing 
the  wave  formation  of  a  linear  and  a  non-linear  computer 
simulation.  Wave  contour  pr  inter-i)lot  s  of  l'ij;ures  9  and 
10  show  the  non-linear  and  linear  initial  \vaves,  respec¬ 
tively,  at  four  kilometers  from  tlie  explosions.  hach  line 
of  tlie  contour  plots  indicates  a  change  in  height  of  two 
meters.  The  trouglis  or  areas  of  the  water  surface  below 
the  normal  free  surface  are  not  represented  liy  contour 
lines,  but  left  blank.  I'he  plots  shov,  how  tlie  non-linear 
effects  produce  a  slightly  smoother  aiul  larger  wave  in  the 
same  distance.  However,  the  linear  simulation  jirovides  a 
very  good  approximation  for  tlie  initial  conditions  of  this 
study . 
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iMon  - 1,  i  iie;i  r  Contour  I’riiiter -Plot 


AloiU’  \i  i;ty  of  I  xpiosions 
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I  1 II  ri.'  1  0  . 
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An  ariay  of  s  iiiiu  1 1  ancons  uiiLlcrwalcr  explosions  is  very 
effective  in  generat  ini^  a  plane  tidal  wave.  i'his  conclusion 
is  based  on  slialiow  water  theor>  with  tlie  explosion  represented 
by  a  raised  cone  of  water  which  produces  waves  with  a  wavelength 
of  about  IhUO  meters. 

The  sensitivity  of  the  wave  generation  process  was  deter¬ 
mined  by  the  height  of  the  initial  wave.  The  wave  heiglit  was 
not  very  sensitive  to  nu  sa  1  i  gnnien  t  s  less  than  .100  meters,  time 
delavs  less  than  5  seconds,  bottom  slope  less  than  0.0004,  and 
ratios  of  yields  in  tlie  range'  of  0.5  to  1.0.  i'lic  wave  height 
sensitivity  to  spacings  of  1.0  ±  .2  kilometers  is  aoprox imat e  1  y 
-11.25  meters  per  kilometers  sjiacing.  Die  yield  sensitivity, 
whicli  depends  on  the  derived  initial  condition,  is  4.55  meters 
per  cube  root  of  the  yield  in  kilotons.  Die  non-linear  effects 
increase  the  size  of  the  tidal  wave. 
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V'l  .  RecoiiiiMC'iKlat  ion 

More  researeli  on  the  initial  eoiulition  is  recoiiiiiiended . 

The  waves  produced  by  the  raised  cone  of  water  have  a  wave- 
length  about  three  times  too  lony  and  tlie  initial  wave  is 
not  preceded  by  a  trouyh  as  normally  expectotl  iUt,'!'  5;  275). 

In  addition,  the  cone  does  not  conserve  tlie  volume  of  water. 

A  new  initial  condition  should  be  ileveloped  which  conserves 
the  water  volume  and  iiroduces  waves  witli  a  shorter  wave  1  en_y  t  li . 
It  is  recommended  that  a  cvliiulrical  cavity  w  i  t  It  its  dejith 
equal  to  its  radius  (  1  (Ml  meters)  and  with  a  lip  of  water 
surrounding  it  be  used  as  a  new  initial  condition. 
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20  abstract  f CoriMriifO  on  r»veff«  ani^  It  ricN  anil  hlantHy  by  M«-*  fc  nimN  rr) 

This  report  considers  an  array  of  simultaneous  underwater 
nuclear  explosions  generating  plane  tidal  waves  on  the  continental 
shelf.  Computer  simulations  are  used  to  studv  this  wave  genera¬ 
tion  process.  They  are  liased  on  non-linear  shallow  water  theory 
with  an  initial  condition  of  a  stationary  raised  cone  of  water. 

The  sensitivity  of  the  wave  height  to  different  I’ields,  eombina- 
tion  of  different  yields,  spacing,  alignment,  timing,  and  bottom  ^ 
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slope  is  del  e  nn  i  iiei-i .  The  lererence  values  used  (or  this  study 
are  a  yield  of  five  kilotons  INi'l  ,  oiu-  kiloiiietei'  spaeiny,  water 
depth  of  100  meters,  aiul  no  hottrnn  slope.  in  aildition  to  the 
sensitivities,  tiie  amount  of  non-lineai  effect  is  illustrated 
by  comparison  of  wave  contour  p  r  i  n  t  e  r  -  p  1  o  t  s  ye/ierated  by  linear 
and  non-linear  shallow  water  tlieory.  The  computer  results 
indicate  the  wave  yeneratiirn  is  relat  ively  ijiseiisitivc  to  real¬ 
istic  variations  in  tlie  |)a  ramet  e  r  s  . 
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